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Abstract 

The formation probability is shown to increase by a few orders of mag- 
nitude if microscopic transport coefficients are used rather than those of 
the common macroscopic pictures. Quantum effects in collective dynam- 
ics are taken into account through the fluctuating force, as exhibited in 
diffusion coefficients for a Gaussian process. In the range of temperatures 
considered here, they turn out to be of lesser importance. 
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1 Introduction 



The description of fusion of heavy ions is commonly performed by studying sep- 
arately the entrance phase and the formation of the compound nucleus 0-0- 
Typically, for the latter process the system has to overcome a barrier and in this 
sense resembles nuclear fission. Both processes may be described by transport 
models involving Fokker-Planck or Langevin equations. One essential difference 
of both cases is found in different initial and boundary conditions. Whereas fis- 
sion can be pictured as the decay of a meta-stable state the formation process 
starts out from some intermediate distribution in collective phase space. Often 
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the latter is simply assumed to be one which moves towards the inner barrier, 
a picture which we will adopt in this paper. Our main aim will be to examine 
the importance of using microscopic transport coefficients which in one way or 
other are governed by quantum effects. To this end we are going to compare such 
calculations with results obtained within the commonly assumed macroscopic 
picture. For thermal fission the influence of microscopic transport coefficient has 
been shown to modify the decay rate considerably [3 |S] . 



2 The formation probability 

We assume the formation of the compound nucleus to be described as a Gaussian 
Markov process over an inner barrier, see Fig.0 In order to allow for a largely an- 
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Figure 1: Sketch of the scenario of this paper. 



alytic approach the barrier will be supposed to be given by an inverted oscillator. 
Restricting ourselves to Gaussian distributions in phase space the probability to 
overcome this barrier is simply given by [Oj 



n(t) 



1 + erf 



erf 



2E„{t) 



(1) 



The q(t) = (q) t specifies the center of the Gaussian, which moves from left to 
right and which within this approximation represents the average trajectory. The 
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^qqif) = ((<? — (l)t) 2 )t defines the fluctuation of the coordinate. Analogous def- 
initions are used for the kinetic momentum P, both for the average trajectory 
P(t) as well as for the fluctuations £ pp (£) and H qp (t). Analytical expressions for 
q{t) and X g(? (t), which may be derived from a Fokker-Planck equation [TUj . are 
given in appendix El Fortunately, when time progresses past a certain "relax- 
ation time" (t ^> 1/|2:~|) the probability (JTJ) approaches a finite stationary value: 

Hoe = lim II(t) = -erfl ^= ) . (2) 

The time dependence of both the coordinate as well as of the fluctuation is 
determined by the frequencies (see appendix EJ 



z ± = — I ±Jl + 4 - 1 ) . (3) 

The Tk in measures the relaxation of the momentum to the thermal Maxwell dis- 
tribution and the (dimensionless) parameter rjf, indicates whether the process is 
under-damped (rjb < 1) or over-damped (rjb > 1). In terms of the transport co- 
efficients for average motion, inertia M, friction 7 and stiffness C, which here is 
negative C = —\C\, these quantities are given by 




M 7 1 Fr~a fA . 

Tkin= v and r]b= ^m = ^^- () 

For strong damping rjb 3> 1, viz r co n 3> Tki n , we have 

z + ^ — = H and >-J_. (5) 

Tcoll 7 r kin 

In this limit only processes are relevant (see (|21|) and (|23|)) which take place on 
the time scale r co n. Later on we will also need the frequency zuf, at the barrier, 



*b=\l l 3 = -F=^=- ( 6 ) 



The dynamics of fluctuations is not only controlled by the transport coefficients 
of average motion but by diffusion coefficients as well. In the classical limit (for 
collective motion) there is only the one given by the famous Einstein relation 
D pp = 7T. In the quantum case and for larger damping there may be another 
one, D qp . Both are to be calculated from the quantal fluctuation-dissipation 
theorem 1 , if the latter only is generalized to the barrier problem by a suitable 

■"■It is for this reason that there are no more than 2 diffusion coefficients 
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analytic continuation [TTJ [TU] . We will address this question in more detail later 
in sect. 



Finally, we like to note that, different to the original papers, we want to use di- 
mensionless coordinates and the conjugate momenta of dimension H. This means 
to apply the following canonical replacement 



ntOb V Mxub 



which for the second moments implies 



Sgq — > T, qq , Spp — -> Spp , S gp — > Sq p . (8) 

nvjh M Wh 



3 Transport coefficients of average motion and 
diffusion 

To distinguish between the coefficients of average motion, M, 7 and C, and 
the diffusion coefficients D^, the notion transport coefficients will henceforth 
be reserved for the former. Later on we are going to compare results obtained 
for microscopic and macroscopic pictures of the dynamics of nucleons inside the 
nucleus. In the first case this refers to quantal motion of nucleons in a mean field. 
Within the Gaussian approach quantum effects in collective motion only show up 
in the diffusion coefficients and, hence, only in the dynamics of the fluctuations. 
With respect to this property the notion macroscopic limit stands synonymous 
for the high temperature limit for which the Einstein relation applies. 



3.1 Transport coefficients of average motion 

Microscopic model: Concerning the microscopic transport coefficients we use 
the suggestion made in jHj 

km ~ 0.6T 2 MeV 1 } 

0.6T 2 n__ 

TcoU " 1 + ^VcLcro MeV (10) 

(with temperature being measured in MeV). These simple expressions have been 
developed to resemble the T-dependence of the truly microscopic values of Tki n (Q, T) 
and r co u(Q,T) obtained by application of linear response theory. Using Q and 
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(fTUj) we calculate the effective damping strength rj b (T) from (HJ) and the frequency 
Wb(T) from ©. As may be seen from Fig.|2]the latter may fairly well be approx- 
imated by tJb(T) ~ 1 MeV/h. In f|T0|) the c macro is introduced as a parameter to 
fix the high temperature limit of friction. In more recent times this is claimed 
to be given not by the wall formula 7 wa n but by about one half of this value, 
see e.g. [H] and further references given there; we adopt this suggestion. It may 
be mentioned in passing that the estimates given in Q and (fTUj) could be ex- 
tended to include the influence of pairing jHj. The latter diminishes friction at 
lower temperatures T < 0.5 MeV but its influence is negligible in the range of 
temperatures used in the present application. 
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Figure 2: Temperature dependence of relevant ratios of transport coefficients. 
Critical damping 7#, = 1 is marked by a dotted line in the bottom right panel. 
Crosses refer to the microscopic computations of [Hj (with pairing included). The 
solid line represents the fits (JHI) and (fTUj) to these data. Macroscopic transport 
coefficients are given by the dashed line [T2] . 



Macroscopic model: In jlj the friction coefficient is taken from the wall-and- 
window formula. Like in jHj we will use 7 = 7 wa ii/2 here. The macroscopic value 
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of the frequency zUb(T) can be calculated from the liquid drop values of stiffness 
Cldm(^) and inertia M L dm; we followed the procedure described in |T2j. The 
ratios Tkm(T) and rjb(T) are then fixed by (J1J). 

In Fig. El we show a comparison of the temperature dependence of the mi- 
croscopic and macroscopic transport coefficients used in this paper. It is most 
important to realize that in the microscopic case collective motion changes from 
under-damped to over-damped at a temperature T « 2 MeV, whereas in the 
macroscopic limit collective motion is over-damped (r^ > 1) even at very small 
temperatures. Please have in mind that with increasing r/& the inertia M becomes 
less and less important, to drop out entirely for the truly over-damped case. For 
this reason the estimates given above should be understood to represent those 
quantities which involve M only for smaller temperatures, up to T « 2 ... 3 MeV 
at best. 



3.2 Diffusion coefficients 

As mentioned already before, the diffusion coefficients are to be calculated from 
the fluctuation dissipation theorem. The latter specifies equilibrium fluctuations 
E^, to which the dynamical ones turn to at large times (see appendix [B] or 
[10 J. To evaluate the for a barrier, as represented by an inverted oscillator, 
requires a suitable analytic continuation involving the change of the stiffness from 
|C| to C — — |C|; details of the procedure may be found in [TTJ [THj - Accounting 
for our presently used transformation (jHj) the relations to the diffusion coefficients 
become 

D m = ^ e q f = > d» = h ™» - ^vt/V) > d pp = (ii) 

Please note that the E^' 6 of the inverted oscillator is negative. 

The application of (fTTJ) becomes particularly simple in the high temperature 
limit. The classic equipartition theorem of Statistical Mechanics implies that one 
has 

\^f\ = Z;f/h 2 = j?-, £^ = for T (12) 

This is easily recognized considering (JBJ) again. As a consequence, the diffusion 
coefficients turn to the Einstein relations, which in our case read 

1 HT 

Tkin ZUb 

For very small damping formulas (|T2|) and (|T3|) keep their form even in the quan- 
tum case if only T is replaced by the effective temperature 

T , = &p cot /k^ (14) 



V£-= = 2S»T, D b J=D'- T =0. (13) 
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see [T3| or JU]- As this requires rjj, <C 1 this limit sometimes is called the zero 
damping limit. It is seen that friction neither appears in the classical or high 
temperature limit (|T2|) nor, by definition, in the limit |T4j) of zero damping. 

At lower temperature this is no longer true for the general situation where the 
size of damping plays an essential role on the temperature dependence of the 
and hence on the diffusion coefficients. This has already been demonstrated in 
ITU] . In these calculations the transport coefficients of average motion simply 
appeared as given, constant parameters. As discussed above, they themselves may 
vary sensitively with T. For this reason the calculations of pTJ have been repeated 
using now the specific transport coefficients for average motion described above. 
The results for and S^' 6 are shown in Fig. El It must be said that the 

calculation of the momentum fluctuations requires regularization of a frequency 
integral, which comes in through the fluctuation dissipation theorem (j2HJ). Also 
with respect to this problem the same procedure has been applied as in [TT] using 
a Drude regularization with uj-q ~ lOro as the relevant cut-off parameter. As 
we have Wb(T) pa 1 MeV/h (see bottom left panel of Fig. |2J) we took the cut-off 
parameter ud = 10 MeV/h after convincing ourselves that the results are not 
extremely sensitive to this choice. The variation with uo^ is more pronounced 
in the case of macroscopic transport coefficients where for uo^ = 5 . . . 20 MeV/H 
the momentum fluctuations differ from our choice by an additional factor 

VP J 

of 0.8 . . . 1.3 at T = 1 MeV. In the microscopic case the corresponding factor is 
only 0.8 . . . 1.1. Paying attention to the fact that the effective damping is quite 
different in both cases, this is in good agreement with the findings of |llj . 

In Fig. |3]the T-dependence of S^' 6 is exhibited both for the microscopic trans- 
port coefficients (left panel) as well as for those in the macroscopic limit (right 
panel), as described in Sect. 13.11 In each case the fully quantum result is com- 
pared with the limits of high temperature (classical) and zero damping. It is seen 
that at very small temperatures the microscopic results are very well represented 
by the limit of weak damping. For Tf^ b the deviations grow faster, a feature 
which might by related to the fact that the regularization procedure was taken 
over from ^JJ without modification. This problem has not been studied any 
further, simply because the effect of the Spp 6 on the transmission probability is 
weakened by other effects, in particular at larger damping, as will be explained 
shortly. In any case, the effect seen at small T simply results from the fact that 
the microscopic damping strength rjb becomes very small (see lower right part 
of Fig. EJ. Let us note that at large temperatures the microscopic and macro- 
scopic diffusion coefficients approach each other, if normalized to hwb(T), simply 
because this is so for the transport coefficients of average motion discussed in 

Sect, mm 
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Figure 3: T-Dependence of the equilibrium fluctuations of coordinate and mo- 
mentum calculated in various limits for microscopic (left) and macroscopic trans- 
port coefficients (right). The inserts show the behavior at higher temperatures. 



4 Numerical results 



In the last two subsections it was shown that a proper microscopic treatment of 
the quantum heat bath has two major effects: (i) the effective damping strength 
r)b becomes much smaller than that of macroscopic models and (ii) different to 
the Einstein relation a second diffusion coefficient D qp shows up. After specifying 
the initial conditions for our numerical calculations we will examine these effects 
separately. The effect of quantum diffusion on the formation probabilities (0) 
and (0) is contrasted with the limits of high temperature (vanishing D qp ) and 
zero damping. This will be done both for microscopic as well as macroscopic 
transport coefficients of average motion. 



4.1 Initial conditions 

In terms of the barrier height and the initial kinetic energy, expressed by 

E b = -V = ^g 2 and = (15) 

through the initial values of the average coordinate go and the momentum Pq, 
the argument of the error function in (j2J) turns into 

goo a/ Ef n / hw b + y/E b / hw b (z~ I w h ) 

\J&M - Yflh/W + {z-/w b y {EM - Y> e q f) - 2(z-/w b ) Y>M/h 

(16) 
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see appendix [X] A very similar result has been obtained by application of 
Langevin dynamics to the inverted oscillator [T4*] . 

For the barrier height specified in (JT5J) we choose E b = 10 MeV. This is in fair 
agreement with the typical heights of the inner barriers used in The initial 
kinetic energy E$ m and the temperature To depend on the strength of friction 
in the approach phase. As in the present paper this problem is not addressed 
in detail we will simply assume the incident energy E over the barrier top to be 
strongly dissipated such that a scenario holds as outlined in Fig. [TJ As can be 
seen there, the sum of the initial energy E and the barrier height E b are split 
into kinetic energy and thermal excitation according to 

£ kin = R(E + E b ) and E* = aT 2 = (1 - R) (E + E b ), (17) 

with the ratio < R < 1 to be fixed later. As an estimate for the level density 
parameter we use the value a ~ A/10 MeV where A is the nucleon number of 
the compound nucleus. In order to keep the computation simple, we assume 
that temperature stays constant at its initial value T = T throughout the whole 
formation process. In Fig. 0] we show the dependence of To on the energy E for a 
system of nucleon number A = 224 and two different assumptions for R. In [TK] 
temperatures in the range T = 0.68 . . . 1.24 MeV have been studied for a heavier 
system with A = 298 nucleons. On the basis of the surface friction model it has 




- R = 0.25 - 

°0 50 100 150 200 

E in MeV 

Figure 4: Dependence of the initial temperature on the incident energy. 

been argued that almost no kinetic energy is left after the approach phase, viz at 
the starting point go ( see f° r instance Fig. 3 of [3] and Fig. lb of 4J). In contrast 
the authors of Ej have suggested damping in the approach phase to be weaker 
than that of surface friction. In this work we do not want to elaborate on the 
approach phase but simply incorporate this uncertainty by comparing the case of 
vanishing initial kinetic energy E^ m = 0, for which the parameter of (|17|) becomes 
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R = 0, with the one of a weaker friction force at large elongations, chosen such 
that R = 0.25. 

To specify the Gaussian of the initial phase space distribution we need the 
widths in q and P. They shall be parameterized in terms of the E^' 6 by intro- 
ducing two positive numbers w and W, 

E pp (t )=wY%*, Z gg (t )=W— E 9P (t )=0. (18) 

The w fully determines the momentum fluctuation and the W the one in the 
coordinate provided the former is given. In case of a strong friction force in the 
approach phase the w will be close to unity, as then the momentum distribution 
will be close to Maxwell distribution ^0] • The W can be said to define a kind of 
measure for the overall width of the phase space distribution. Indeed, requiring 
W > 1 warrants the uncertainty principle 

((?- (l)t ) 2 )t ((P ~ (P)t ) 2 )t a = X qq (t ) E pp (t ) > h 2 /A (19) 

to be fulfilled initially. In the limiting case W = 1 the system starts out of the 
most narrow distribution still compatible with quantum mechanics. As we want 
to examine the importance of genuine quantum effects vanishing initial fluctua- 
tions £ M j,(t ) = are prohibited for our studies. One must have in mind, however, 
that in reality finite fluctuations will have been built up in the approach phase, 
which may not always be represented by simple Gaussians. In our numerical cal- 
culations, the which enters the initial conditions (|18j) is always taken to be 
the quantum value (see Sect. 13. 2|) . even for studies of the time evolution within 
high temperature limit. 

From the results displayed in Fig. El it can be seen that the E^; 6 of the zero 
damping and high temperature limits may differ considerable from those of the 
full, quantum results. This difference is particularly large for T3> b when macro- 
scopic transport coefficients of average motion are used. One might be inclined 
to believe that at small temperatures this might have a large impact on the for- 
mation probabilities H(t) and Iloo. However, one must be aware that not only the 
E^; 6 contribute to the E^ of the denominator of (|TH|l but the initial fluctuations 
£71!/ (to) as we ll- To reveal these effects we address in Fig. El the T-dependence of 
the E^ . Shown there is the ratio of its values in either the high temperature limit 
or the zero damping limit to that for the full calculation, for the choice w = 1 and 
W = 1 — for other cases this effect is even smaller. It is seen that the values E^ 
obtained in these calculations are quite close to each other. Note, please, that in 
our studies only temperatures T > 0.6 MeV are accessible (see Fig. Hj). As only 
the square root of the E^ appears in (J2J) one may already deduce here that the 
three calculations will lead to results which do not differ too much. 
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Figure 5: Temperature dependence of the for the zero damping and high 
temperature limit, both divided by the general result: left panel for microscopic 
input, right panel for macroscopic input. 



4.2 Formation probabilities for microscopic input 

We plot the time evolution of the formation probability using the microscopic 
set of transport coefficients and the fully quantum calculation of the diffusion 
coefficients in Fig. |U1 Results are shown for an incident energy above the barrier 
of E = 10 MeV and for vanishing (R = 0, Eq 1u = 0) and finite initial kinetic 
energy (R = 0.25, Eq 1u = 5 MeV). For the parameter w, which according to 
(fT%|) specifies the initial width in the momentum, three different values are used: 
w = 1/2 (more narrow than thermal equilibrium), w = 1 (thermal equilibrium) 
and w = 2 (wider than thermal equilibrium). The initial distribution in the 
coordinate is specified by the two values W = 1 and W = 3. In most cases the 
initial space distribution is almost totally located on the left of the barrier. For 
the choice w — 1/2 and W = 3 a small part of the distribution is on its right 
already at the beginning. In this case the formation probability II(i) does not 
start from exactly zero. 

In Fig.Elwe show examples corresponding to small effective damping strengths: 
7] b = 0.26 and r] b = 0.18 at T = 0.94 MeV (R = 0) and T = 0.82 MeV (R = 0.25), 
respectively. As this corresponds to under-damped motion the initial kinetic 
energy plays an important role on the formation probability. Its stationary value 
IIoo increases by two orders of magnitude in going from R — to R = 0.25. 
Repeating the same calculations at larger energy E = 50 MeV reveals additional 
effects. Enlarging R from to 0.25, which is to say for increasing initial kinetic 
energy and slightly decreasing temperature (Ey m = 0, To = 1.6 MeV in the case 
R = and Eg* = 15 MeV, T = 1.4 MeV in the case R = 0.25), the asymptotic 
value IIoo of H(t) is reached quicker, independently of w and W. The influence 
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of the parameters W and w on the value of Iloo is less important than the one 
of R: Varying either w or W mostly implies a change by less than one order of 
magnitude. 
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Figure 6: Time evolution of the formation probability H(t) for microscopic trans- 
port coefficients and quantum diffusion coefficients, for different choices of the 
parameters specified in the inserts; the incident energy is E = 10 MeV. Note 
that the scaling in the upper and lower panels differs by two orders of magnitude. 
For more details see text. 



Next we discuss the energy dependence of the asymptotic value Iloo. This is 
done best by introducing an "effective barrier height" by the following argument: 
Within the Gaussian approximation the transmission across the barrier is largely 
governed by the average trajectory. As can be seen from (J2J) Iloo becomes equal 
to 1/2 when the trajectory just hits the top of the barrier at = 0. By 
inspecting the numerator of (|16jl . and recalling from eq. ([5|) the value of z~ , one 
easily recognizes that the trajectories may be classified by their initial kinetic 
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energy Eq 1u in comparison with the "effective barrier height" (see also jHElEl) 



B, 



eff 



l + r] 2 b +r] b \ E b > E, 



-6 • 



(20) 



Trajectories with Eq 111 < B c g are simply reflected implying the Iloo to be smaller 
than 1/2; only those trajectories with Eq 1u > B c g are able to cross the barrier 
and the Iloo exceeds 1/2. Fig. shows the functional dependence of Iloo on 
the initial kinetic energy scaled to the effective barrier height, for R = 0.25, 
w = 1 and W = 1 (left) and W = 3 (right). For several reasons the incident 
energy is restricted to the interval E = . . . 200 MeV, with a corresponding 
temperature of To = 0.6 . . . 2.6 MeV. First of all, at larger temperatures the 
survival probability becomes too small. Secondly, as mentioned already before, 
there the motion becomes strongly over-damped, such that formulas which involve 
the inertia become less and less trustworthy. At small incident energies E (small 
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Figure 7: Formation probability Iloo using microscopic transport coefficients as 
a function of the ratio of the initial kinetic energy (fTTj) to the effective barrier 
height fl2()D . For we compare the results applying a quantum calculation 

(comp. pi]) and the limits of high temperature ((12)) and zero damping f!12D with 
(I14J) . Note that the incident energy changes along the curve: E — . . . 200 MeV. 
The value E = 10 MeV corresponds to Eq 1u / B eS = 0.35; for details see text. 



T ), for Iloo a choice of W = 3 implies a factor of the order 5 over the case of the 
most narrow initial distribution which still is compatible with quantum mechanics 
(W = 1). This increase of Iloo with W rapidly diminishes with growing E. In 
addition to the full quantum calculation of the diffusion coefficients we show in 
Fig. [7| the results of the limits of zero damping and high temperature. Both 
approximations are of about the same quality. For the chosen example quantum 
effects are seen to be small and become largest at intermediate temperatures 
T = 0.7 ... 2 MeV (E = 5 . . . 120 MeV). Both observations confirm the conjecture 
made already in connection with Fig. 
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As a remarkable feature we observe in Fig. [7| a back-bending of the curves 
representing at E^ m / B c r pa 0.5. It can be explained as follows: The initial 
kinetic energy depends on E linearly (see (fTTjl ). while the independence of the 
effective barrier height (ffllj) is more complicated and comes in via rjb{To(E)). As 
a function of E for R = 0.25 the ratio E^ m / B e Q develops a maximum and a 
minimum (the latter corresponding to T ~ 4 MeV; which is beyond the region 
considered here). At energies far beyond this range IIoo would approach 1. The 
final formation probability also exhibits a shallow maximum near the turning 
point. This behavior might have some influence on the reasonable choice of 
incident energies in an experimental setup for heavy ion collisions, because the 
survival probability decreases with increasing excitation energy E* and therefore 
smaller E is more favorable for survival. 



4.3 Formation probabilities for macroscopic coefficients of 
average motion 

In Fig.|Hlwe show IIoo as a function of the ratio E$ m /-B e ff on a logarithmic scale for 
E = . . . 200 MeV and R = 0.25, w = 1 and W = 3 using macroscopic transport 
coefficients now. The back-bending found in Fig. is absent here. Quantum 
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Figure 8: Same quantities as in the right panel of Fig.0for macroscopic transport 
coefficients. The value E = 10 MeV corresponds to EQ m /B cfi = 0.0036. Note the 
logarithmic scale. 

effects of the collective motion qualitatively behave as discussed in Sect. 14.21 
Quantitatively at T = 0.8 . . . 1.0 MeV they are even smaller than in the case of 
microscopic transport coefficients, comp. Fig. 0D The small differences between 
the full quantum calculation (of the diffusion coefficients) and the limits of zero 
damping and high temperature are hardly visible on logarithmic scale. 
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Due to the fact that for macroscopic transport coefficients collective motion 
is strongly over-damped (jjf, = 5.9 . . . 4.9 for T = 0.6 ... 3 MeV) the effective 
barrier height (|20)l is very large now: 5 e fr = 1.4 GeV. . .980 MeV. Mind that 
for microscopic transport coefficients the same quantity is in the range B e g = 
15 . . . 180 MeV only (comp. E b = 10 MeV). This leads to ratios Ef n /B eS of 
about one order of magnitude smaller than for microscopic transport coefficients. 
This ratio crucially influences the size of the formation probability Iloo. Conse- 
quently in the case of macroscopic transport coefficients Iloo is several orders of 
magnitudes smaller than for microscopic ones. This effect increases with decreas- 
ing E. Compared to this loss of orders of magnitude in Iloo due to the replacement 
of microscopic by macroscopic transport coefficients, quantum effects in the col- 
lective motion due to the type of approximation applied in the calculation of the 
Tfff are of minor importance only. 
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Figure 9: Same quantities as in Fig.EJfor E — 10 MeV and macroscopic transport 
coefficients. 

The time dependent behavior of the formation probability Tl(t) is shown in 
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Fig. El for the same initial conditions as for Fig. |U1 As mentioned previously, 
now the effective damping of the collective motion is clearly in the over-damped 
regime: rj b » 5.8 both at T = 0.94 MeV (R = 0) and T = 0.82 MeV (R = 
0.25). The time scale for the rise of IT(£) is larger by a factor of the order of 10, 
as compared to the microscopic case and completely independent of the initial 
kinetic energy Eq 1ti , viz the initial momentum P . We have checked that different 
to the case of microscopic transport coefficients this is still true at E = 50 MeV. 
This is in accord with the fact that for rjb ^> 1 one is in the Smoluchowski limit, 
for the momentum drops out entirely. For the same reason also the stationary 
value IIoo is much less sensitive to the initial kinetic energy than in the case 
of microscopic transport coefficients. Whereas in the microscopic case at E = 
10 MeV we gained two orders of magnitude going from R = to R = 0.25 the 
increase of now is invisibly small. Using the large energy E = 50 MeV the 
corresponding increase is at most by a factor of 2 whereas in the microscopic case 
we gain two orders of magnitude, too. 



5 Summary and conclusion 

Concentrating on the formation probability we have been able to demonstrate 
the importance of microscopic transport coefficients of average motion. As com- 
pared to those of macroscopic pictures like that of the liquid drop model and wall 
friction, they may increase the formation probabilities by a few orders of magni- 
tude. This is especially important at smaller temperatures where the dynamics of 
nucleons in a mean field is dominated by quantum effects and where "collisional 
damping" is suppressed [TU] . For collective motion this implies a moderate to 
weak damping. 

In contrast to this, quantum effects in collective motion would become relevant 
only at even smaller temperatures. In a sense, this is in accord with the findings 
in [7j for nuclear fission. Since at small T friction is suppressed quantal diffusion 
coefficients may be estimated within the so called "zero friction limit". This 
feature applies only for the microscopic picture of nucleonic dynamics. In fact, in 
this paper it was the first time that microscopic transport coefficients have been 
used when calculating the diffusion coefficients for collective dynamics through 
the quantal fluctuation dissipation theorem. Because of the delicate temperature 
dependence of the transport coefficients this procedure is vital to get decent 
answers. 

Another issue of our studies was to examine the relevance of the initial condi- 
tions for the phase space distribution. Whenever the whole process of fusion is 
split into parts where the entrance phase is treated differently from the formation 
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of the compound nucleus uncertainties for the initial conditions of the formation 
process are more or less inherent. Evidently, such a situation would only be given 
if the intermediate step would correspond to a kind of quasi-equilibrium; in re- 
mote sence in analogy to N. Bohr's hypothesis of a two step process. Indeed, 
as we have seen, for a more quantitative calculation of the reaction cross section 
it will be inevitable to follow the whole process, entrance phase plus formation. 
It is only then that the "initial conditions" for the latter can be known better. 
Unfortunately, however, it is only fair to say that a decent microscopic picture 
for the first stage is still lacking. 



Acknowledgments: The authors would like to thank F.A. Ivanyuk for critical 
and fruitful discussions and D. Boilley for helpful comments. 



A Average coordinate and fluctuations 

Using dimensionless coordinates and the conjugate momenta (J7J) the solution of 
the equations of motion for the average trajectory reads in the case of the inverted 
oscillator [Til] : 

= 7T?r ( w bPo/h- q z~) exp(z + (t -t )) 
+ 77?r(-^bPo/fi + q z + ) exp(z~ (t - t )) (21) 
Here go = <?(^o) and P = P(to) are the initial coordinate and momentum respec- 



tively and the frequencies z ± are defined in 0; £& = (2rki n ) 1 \J 1 + rj b 2 . For long 
times t 3> l/k _ | (EH) reduces to 



?(*) ~ i^bPo/h - qoz-) exp(z + (t - t )) . (22) 



In dimensionless coordinates the fluctuations of the coordinate of the inverted 
oscillator read (comp. [TU]): 

2 2 

£ OT (t) = M^p e*v{2z-{t-to)) + A 2 ^ exp(2z + (t-t )) 

+A 3 eM(z- + z+)(t - to)) + Y> e q f (23) 
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Here the abbreviations 
A x = 

A 2 
A 3 

and 



z 



[z^ - z' 
z + 



(z+ 


-z-y 


z + 


+ z~ 


(z+ 


- z-y 






(:! 











z+ 






2 




c + 


- 




w b 




\zu b 




z~ 






2 




c + 


[- 




W b 




\w b 



z + z^ z z 
c + 



w b 



b + a 

b + a 

f 

-b + a 



(24) 

(25) 
(26) 



(E OT (t ) - )//» \ 

(E gq (t )-Z% b )h (27) 

Sgp(to) / 

have been applied (note the small differences in the definitions compared to ^0] 
that make the expressions more symmetric). The long time limit of (j2Hj) reads: 



— 



E M (t) ~ A 2 jr± exp(2z + (t-t )) 



(28) 



B Implications of the fluctuation dissipation the- 
orem 

For a stable system the equilibrium fluctuations of some quantity A^ can be 
calculated from the quantal fluctuation dissipation theorem in its normal form: 

^ = hj — coth — x%(u) (29) 

Here, x'L * s the dissipative part of the response function that describes the vari- 
ation of the average of A^ to a small change of the "external field" a v which 
couples to the system through an interaction 5H = A v a v (t). For the present 
case, the average (A^) either refers to the coordinate q or the momentum p. 

Unstable modes, which appear in the barrier region, may be treated by analytic 
continuation [H], changing a positive stiffness into a negative one. In this way 
the stationary fluctuations of the damped inverted oscillator (of frequency w b 
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and effective damping r] b ) become: 



gq 



S- 6 = ^ 2 (|S^| + ^^x 



KT kin ZU b 

x (uU) - cjO'+ 1 ))(cjO') - cjO'+ 2 )) ' + 1 ~^rrT I I ( ' J, L) 



Eg' 6 = (32) 

Here iuq ^> U7& is the so called Drude cut-off frequency necessary to regularize 
the cj-integral in (J2HJ) in the case of E^-' 6 . We have put zu D = 10 MeV/ft in this 
paper and the frequencies u>^> are the solutions of the equation 

-ml -iul b - + - a; 2 + i = (33) 

V W D T kin J W D 

(periodicity is assumed: u^ +3 ' = uj^). 
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